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SUMMARY 

An exact solution of the problem of designing an airfoil with 
a prescribed velocity distribution on the auction surface in a given 
uniform flow of an incompressible perfect fluid is obtained by 
replacing the boundary of the airfoil by vortices. By this device, 
a method of solution is developed, which is applicable both to 
isolated airfoils and to aiirfoils in cascade. The confoimal 
transformation of the designed airfoil into a circle can then be 
obtained and the velocity distribution at any angle of attack 
computed. Numerical illustrations of the method are given for the 
airfoil in cascade. 


INTRODUCTION 

The problem of increasing the output per stage in axial-flow 
compressors and turbines involves the use of high -solidity (closely 
spaced blades) stages of highly cambered blades. In addition, the 
velocity distribution must be carefully selected as a function of 
arc length along the airfoil (blade section) boundary in order to 
avoid flow separation or excessively high local velocities. 

Several methods are available of obtaining an airfoil with a 
prescribed velocity distribution when placed in a fixed position in 
a uniiorm stream. The methods that lead to theoretically exact 
results are based on conformal mapping theory. (See references 1 
and 2.) In reference 3, Mutterperl extended the method of conformal 
mapping to solve the problem of computing a cascade of airfoils with 
prescribed velocity distribution but, for cascades with closely 
spaced or highly cambered airfoils, this procedure becomes very 
cumbersome. Approximate solutions have been obtained by placing 
singularities such as vortices, sources, and sinlcs in a unifonn 
stream. Cascades of airfoils can also be computed by distributing 
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such singularities periodically throughout the region of the cas- 
cade, as descrihed hy Ackeret (reference 4). Because fnese vor- 
tex methods are not exact, a method has been developed with the 
vortices on the boundaries of the cascade airfoils. This results 
in a theoretically exact solution without the computation difficul- 
ties encountered in conforaal mapping methods for highly cambered 
airfoils or closely spaced cascades. Furthermore, for the same 
accuracy in computing the airfoil shape, this vortex method re- 
quires the computation of fewer points than the method of conformal 
mapping because these xx)ints may be arbitrarily placed on the airfoil. 
The method may be applied to isolated airfoils and to airfoils in 
cascade. For the cascade, the inflow and discharge velocities and 
a velocity distribution on the surface of an airfoil are given and 
the shaxje of the airfoil is determined. In some cases, the spacing 
of the blades is preassigned, which places a condition on the as- 
sumed velocity distribution. Before the airfoil shape has been 
evolved, the velocity distribution may be computed foi’ any angle 
of attack by the method described in appendix A. 


THEORY 

Outline of Method 

In ref'-rence 5, it is demonsti'ated that the two-dimensional 
potential flow about a body in a uniform stream can be represented 
bj’’ substituting for the body a sheet of vortices along its boundary. 
The vortex strength per arc length at any point is equal to the 
magnitude of the velocity at that point. A proof of this relation 
for the case of the cascade is given in appendix 3. The problem 
of finding a shape with a prescribed velocity distribution when placed 
in a stream can then be stated: given a vortex distribution, to find 

a contour which satisfies the condition that it will be a streamline 
in the flow field induced by the uniform flow and the vortices 
distributed on the contour. 

The procedui'e of finding the shape begins by cboosi.ng an 
approximate shape and distributing the vortices on it. The stream 
function of the flow induced by the’ vortices and the uniform stream 
is computed at points on the boundary of the assumed shape. If 
this stream function is constant, the assumed shape is correct. 
Variations of the stream function are a measure of the deviation of 
the assumed shape from the correct one. These variations are used 
to distort the original shape into a new shape whose stream function 
is more nearly constant. The process is repeated until the 
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variations 110001116 negligible. In the process of shape adjustment, 
the velocity is altered on the pressure surface. 


Derivation of Equations for the Stream Function 


Isolated airfoil. - The complex or reflected 

(which is the derivative of the complex potential 

induced at the point z by a vortex of strength 

z is 
o 


velocity w'(z) 
function w(z)) 
k located at 


w'(z) 


k 1 
2ni z-Zq 


A summary of the main sjimbola used in this report is given in 
appendix C. 


The complex velocity, w'(z) induced by a uniform, stream with 
complex velocity w^ ' and a distribution of vortex strength per unit 
length 7 (Zq) along a curve with coordinates Zq, is 


w'(z) = Wy' 


1 r 7(zo)d3o 
^ 2niJ z-Zq 


( 1 ) 


where dSg ia the element of £irc length along the curve. The complex 
potential at the point z is the integral of w'(z) with respect 
to z, namely 

w(z) = zw„' + 7 (zg) log (z-Zq) dsQ (2) 

From reference 1 (notation modified) 


l(Zo) = w'(zo) dzo = dw(zo) = dcp(zo) + id\^(zo) 


where 

/ 

cp velocity potential, E[w(z)] 
y\i stream function, I[w(z)] 

When equation (2) is applied to obtain the complex potential 
function at any point z in the flow field, the integration must be 
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carried out along the bovmdary of the tody. Because this curve is a 
streamline, d\|/ = 0 and, therefore, eq,uation (2) becomes 


w(z) 


ZVv 


+ 



log (z-Zq) d(?(zo) 


(2(a)) 


The imaginary part of eq.uation (2(a)) is the stream function at the 
point z. 


\|/ = - xVy + yV^ 



log y(x-Xg)2 + (y-yo)2 dcp(zo) 


(3) 


where 


Wy ' complex velocity of uniform stream, for isolated airfoil, 

- "y) 

z coordinate of point where stream function is computed, (x + iy) 

Zq coordinate of point where vortex is located, (xq + iyo) 

Vx x-component of uniform stream velocity V 

V y-component of uniform stream velocity V 

«/ 

It is convenient to use the arc length along the airfoil as a 
parameter. If (x, y) is a point on the airfoil boundary, then s 

will denote the arc length there; similarly, s^ will denote the 

arc length at (x^, y^) . The vortex at on the airfoil influences 

the stream function at the point s on the airfoil. The stream 
function induced at (x, y) by a vortex of unit strength at 

(xq^ yo^ 


- ^i(z, Zq) = - 




(4) 


A plot in the (x, y) plane of curves for a constant fi(z, Zq) 
consists of concentric circles with center at (x^, y^) . 

The velocity at the point s^ on the airfoil is the directional 
derivative cp’(Sq) of the potential along the streamline. If the 
velocity along the airfoil has been specified and an airfoil shape 
has been assumed, the resultant stream function along the boundary 
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of the airfoil can he approximated hy using the approximate shape in 
evaluating the integral - 


ilia) 


pi 

■/ ^l(®^ ®o) <=P'(So'^ <^^0 
<J0 


where 

\l/^(3) stream function at (x, y) due to uniform stream, 
(-xVy + yV^) 

I total arc length of airfoil 


(5) 


All variables are expressed in terms of the arc-length parameters s 
and Sq. The integral in equation (5) can be evaluated either 
numerically or graphically over the entire range of integration 
except in the region where a = s-Sq is small, for in this 
region becomes infinite. This portion of the integral ' 

can be evaluated by approximating the airfoil boundary by a line 
segment . Then, 


The prescribed velocity can be given in this region, which may be 
defined by s-a s^ s+a, by a Taylor expansion as a function . 
of about the point s . 


cp'(s^) =cp'(s) +cp"(s) (Sq-s) + (Sj^-s)^ + , . 

where the primes indicate derivatives with respect to s. The 
integral is then 


Pa+a 

'-'s-a 


fl(Si'So) cp'(so) dso 
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In most cases it is necessary to use only tlie first, term in equa- 
tion (6). The same type of approximation can he used to evaluate a 
portion of the integral if the opposite side of the airfoil comes 
in the neighborhood of the point (x, y) . 

A more general equation applicable to a segment that does not 
pass through s is: 


cp'(s„) d(s„) 


= ^C'log (h^+c^) - b log (h^+b^) - 2(c-b) + 2h tan 

+ l^(h^-t-c^) log (h^+c^) - (h^+b^) log (h^+b^) - (c^-b^) 


1 h(c-'b) 
h^+bc 


T. 


L'J l jll log (h^+c^) - b^ log (h^+b^) - I (c^-b^) 
3’ L ' ■ 


+ 2h^ (crb). 2h^ tan 


u3 l.„-l Mc-b) 


h?H-bc 


■ • e 


(6(a)) 


where 'h is the perpendicular distance from s to the segment^ 

3^ = ’.p " locates the foot, of the segment, (p+b) and (p+c) are the 

limits of the integration of Sq, and approximately, 


" ■' ' ' ' v2'...,' 

cp'(p) = cp'(p+b) - b cp' '(p+b):’.+ ,cp.' ' ' ( 


cp"(p) =cp"(p+b) -bcp”'(p+b) 

; cp"-(p) = cp"'(p+b) 

Equation' (6(a)) may be used when the line segment is not of equal 
lengths on either side of the perpendicular foot or when cp'(s) or 
its derivatives are discontinuous at either (p+b) or (p+c). If 
a = c = -b and h = 0, equation (6(a)) reduces to equation (6), 
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The size of a, h, or c is determined hy the requirem.ents that the 
segment in question "be nearly straight (the approximation is of the 
second degree) and that cp'(g^) be accurately represented by a 
Taylor series expansion of few terms. 

Airfoils in cascade. - The expression for the complex potential 
for the flow about a cascade of airfoils is derived in appendix B. 
The notation is defined in figure 1. The formula, which corresponds 
to equation (2(a)) for isolated airfoils, is 


w(z) 


= zv, 


m 




2ni 


r 


log 



d'‘p(Zo) 


(7) 


where 

d distance between successive airfoils in cascade 

mean stream velocity, whic;i is one half the sum, of com.plex 
(reflected) velocities upstream, and downstream, of cascade. 


The mean velocity Wjjj ' corresponds to the uniform velocity w^ ' of 
the isolated airfoil flow. 

The term zWj..^^' is the complex potential function I’esulting 
from the mean flow. In the integral, the element dcp indicates the 
vortex element strength and -log [sin (jt/d) (z-Zq)] represents 
—the— complex-potentlal-at— the— point — z — caused-by— an-inf-ini-te-row— of— 
unit vortices at z„ ± nd where n=0, 1, 2, . . .. The imaginary 
part of equation (7) is the stream function. 


J r>SQ=Z 

f2(s, Sq) dcp(sQ) 

So=0 


(8) 


where 


^2^®^ ®o) 



sin^ I (x-Xq) + sinh^ ^ (j-Jq) 


is expressed in arc-length parameters and '4/jj^(s) is the stream 
function at (x, y) induced by a mean stream whose complex velocity 
is Wjjj'; that is. 
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The values of 



for various values of fg are 


given in table I, A plot of x-Xg and y-y^ for constant values 
of f2(z,, Zq) is shown in figure 2. These curves may be interpreted 
as the streamlines of the flow induced by an infinite row of vortices 
of unit strength located at the points (x^ind, y^), where 

n = 0, 1, 2, . . In the region of a vortex, the streamlines are 
nearly circles; that is, the flow is nearly that induced by an 
isolated vortex. At a distance from the vortex row, the streamlines 
are parallel lines, as in the flow pattern induced by a continuous 
uniform distribution of vorticity along a straight line instead of 
a row of discrete vortices. The velocities on the two sides of such 
a vortex line are of equal magnitude but opposite direction. 


This behavior of f2 for large 
[(y-y^)^ + (x-Xo)2] 


|y-yp| 


and also for small 
x-x.. 


and 


y-yo 


can be described as follows; When both 


are small. 


f 2 (z, Zq) « ~ log ^^(x-Xo)2 + (y-yo)^J 


(9) 


Which differs from z^) only by a constant. For large values 

^ I y-yo I . ^"^o 

of irrespective of ~ — and a constant term, 


f2(z, Zq) « (10) 

which is the stream function of a uniform stream parallel to the 
x-axis . 

Equation (G) can be used for computing the stream function along 
the boundary of an airfoil in cascade just as equation (5) is used 
for the isolated airfoil. The integral over the range in the neighbor 
hood of the point s is obtained by using equation (9) for f 2 (^> So) 
The result, derived in the same manner as equation (6), is 
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It 



3 ^) Cp-(S^) dS^ 



Cp' (s) 



+ a 


3 cp ' ' ' (a) 


3.' 




(6(b)) 


The more general equation (6(a)) is modified for cascades by 
multiplying the argimient of all logarithms by the factor it^/d^. 


Adjustment of Shape 

If the stream function for the assumed airfoil has been computed 
and has been found to var 3 '’, the shape must then be adjusted to give 
a more nearly constant stream function. The shape changes are made 
by rotation of the body^ plus displacement of the individual points 
normal to the mean stream. A succession of such changes results in 
the most general type of shape distortion. The rotation is used to 
place the front stagnation point properly. 

Notation of the airfoil. - In the formula for computing the 
stream, function of an isolated airfoil, the contribution of a vortex 
element at (xq, y©) to the stream, function of a point at (x, y) 
is dependent merely on the distance between the two points. 
Consequently, if the entire airfoil is rotated, the effect of the 
-boundary— vorti-ees-on—t-he-stream-functl-on—at-any-po-int—on-t-he-a-irfoi-1— 
boundary will not change. The effect of the blade rotation on the 
stream function along the boundary is therefore determined by the 
change in relative position of the points in the uniform stream. 

The first adjustment in shape is a rigid rotation of the airfoil in 
order to obtain a more nearly constant stream function along its 
boundary. 


If the airfoil is rotated through an angle P, the stream 
function is so changed that \J/(s) is a function of p and s and 
may be written p) . Wlien p = 0, '^/(s, 0) is the original 

stream function before rotation. After rotation the new stream 
function \l/(s, P) may be expanded in a Taylor series about the 
point p = 0, 


vKs, 


p) = >j/(s, 0) + p 


^^l'(a, P) 
Sp 




P=0 
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Only the first two terms in this series will he used because it is 
assumed that 3 is small. The angle P is to be determined for 
the minimum mean-square deviation of the stream function from its 
mean value. Because the object of the rotation is essentially to 
adjust the shape of the nose, the rotation might also be made to 
reduce the root-mean-square deviation of the stream function to a 
minimum .for a portion of the shape including the nose. 

The mean value of the stream function at any angle P is 

0) ,+ p ^•^(s, P^ ^ ds , (11) 


\i/(3) .= jP ^{a, p) ds = jj 


The difference between the new stream function \|/(s, P) and its moan > 
value \KP) is squared and integrated to obtain a measure of the 
variation of, (j/Cs, p) . from the mean value at the new angle . The 
condition for obtaining a minimum root-mean-square deviation by 
adjusting P is , 




. " Jo 


\l/(s, 0) + p - 'K3) 


2 ' 


^ds 


: 2 \l/( S, 0)- + p \1/ ( P }: 


d-4/(s,' 0) - ^MjPi 
dp "dp 


ds (12) 




d\j/(3, 0) : 


dp- 


4^(s,; 0)' + ,p: 


ds 




^0 ^ 


ds = 0 


(13) 


The second integral vanishe_s by. virtue of equation (ll), which may 
also be used to, eliminate ■^/(P),:., from, the remaining term.. The solution 
for- ,P is- 


P = 


JO 


%(s) 


d\|/(s , 0 ) 

:. ■, . .dP • 


■ ir ri /; . T d\i'(s, 0). , . 


(14) 


1 

I 



i2. pi 

.Jo , . 

1. JO 


d\|/(s . 0) f 

. dP J 


ds 
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In order to apply equation (14), d\l//d!3 must te known at points 

along the "boundary of the airfoil. In the case of the isolated air- 
foil, the contribution of the vortices is unaffected by the rotation 
and therefore, 



If the airfoil is rotated about the point (x_, y„), equa- 
tion (15) becomes 


H = cos p[^ (x-x^.) + (y-y^) 

+ sin p [(x-Xq) Vy - (y-y^) 


(16) 


where (x, y) are the coordinates of the point before rotation. 
For small values of P, equation (16) reduces to 


0 = (x-Xp) + (y-y^) Vy (17) 

The choice of (x^,, y^) will have no effect on the results in this 
case. 

When the ajrloil In cascade is rotated., — .the_change — a.-n— the— posi- 

tion of tjie vortices of. the adjacent blade must be considered. In 
the case of the isolated airfoil, it was unnecessary to consider 
the change in position of the vortices because the influence of a 
vortex (equations (5),. (4), and (5)) depended on the function f^, 
which is constant on circles. The influence of the vortices on 
the airfoil is, therefore independent of direction. Because the 
fg contours are not circles, the rotation in cascade 'does have an 
effect, which is approximated by considering all closed fg curves 
(f 2 < 0 ) as circles in order that the effect of all vortices 
in the region fg<0 may be neglected during rotation. The 
effect of all vortices in the region f 2>0 is estimated by ■ 
assuming that all the fg contours for f 2>0 are straight lines 
uniformly spaced. The flow corresponds to that between two infinite 
straight parallel vortex sheets of uniform strength per unit length. 
This flow induced by the vortices in the region fg>0 is in the 
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x-direction, and the direction of the flow induced hy-the vortices 
for which is opposite in sense to that induced hy the vortices 

for which YQ<y> 

As the point being considered is changed, the regions for 
fg^O, yo>y> ^2^^^ yo<y win include different sections 

of the blades, and hence different vorticity, with the result that 
the x-velocity component induced by the vortex sheets will vary 

with the point under consideration. The algebraic sum of the 
x-component of the uniform flow velocity aind the variable x- velocity v^ 
induced by the vortices in the region f 2>0 is to be used like the 
velocity component in rotation of the isolated airfoil (equa- 
tion (17)), The quantity in equation (17) is replaced by the corre^ 

spending V^. ^ (= + v^) . The vortex strength per unit length at any 

point on the airfoil is equal to cp’(so) and, therefore, from equa- 
tion (10) the x-component of the velocity induced by the vortices is 





where the integration is carried out over the portion 


of the airfoil where f 2 (s, Sq)>0, A distinction must be made between 
the two regions yo<y aii-d yo>y because the induced components have 
opposite direction. 


The computed result of rotating an airfoil in cascade depends upon 
the choice of (xq, j^) because d\l//dp is only approximate. In order • 
to minimize the error involved, values of, dAj/ZdlB are reduced by 
choosing (x^,, y^,) as the centroid of the vortex distribution on the 
airfoil. If the improvement in the mean-square deviation of 4; is 
small compared with its original value, it may be preferable to omit 
the rotation of. the airfoil because of the error inherent in the computa- 
tion. The decision should be made chiefly on how 4/ varies at the air- 
foil nose and whether it is approaching a constant value in this region 
with successive corrections of the shape. 

Distortion of the shape. - The stream function computed after the 
isolated airfoil has been rotated will, in general, still vary along 
the boundary.. This variation can be reduced by distorting the shape of 
the airfoil. If the distortion is small, the change in distance between 
einy two points on the boundary will be small, although the change in the 
direction of a segTaent Joining those points may be considerable. The 
effect of the distortion on the contribution to the stream function of 
the vortices on the boundary. is consequently neglected. The largest 
effect of the distortion will be to change the position of the boundary 
points in the uniform stream. The airfoil is therefore distorted in such 
a manner that the change in the contribution of the unifona stream to the 
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stream function will eliminate the variations in stream function. 

For points directly opposite each other on the airfoil^ the change 
in distance will he of the same order of magnitude as the distor- 
tion. Consequently, distortions that result in change of thickness 
of the airfoil converge very slowly because of the Inaccuracy of the 
fundamental assumption on which the correction is based. 

Thus, when the stream function along the boundary of the 
isolated airfoil is known, some number is arbitrarily chosen as the 
desired constant value of the stream function. If = vj/ - is 
the difference between the computed stream function at a point and 
the desired constant, the point is moved a distance •• A\l//V perpen- 
dicular to the dii’ection of the mean stream, where the direction of 
increasing uniform stream function is taken as positive. The air- 
foil in a cascade is distorted in the same manner, using the varying 

resultant local mean stream velocity corrections 

are made with ‘4/ eq.ua! to the mean value of \1/ on the airfoil. 

COMPUTATIONAL PEOCEDUEE FOE CASCADES 
Choice of Velocity Distribution 

Several factors influence the choice of the velocity distribu- 
tion for which an airfoil is to be found. Especially in rotors, 
sturdy blades are required. Long thin tail sections must be 
avoided and where high rotative speeds and stresses occur, overhang 
of thin sections is likely to induce blade failure. The radial 
"d rs tr ibut f on L>f “a^irf oil” cross - sect ioria.l area is also fundamental in 
determining the blade-root stresses. Overhang can be reduced 
by proper choice of the velocity diagrams for the sections, but the 
other factors are influenced chiefly by the thickness of the section. 

The desired thickness may be attained by first assuming a blade 
shape and spacing and using the stream-filament method of reference 6 
to compute the velocity distribution over a portion of the airfoil 
that determines the thickness. The spacing may be regarded as fixed 
but the curvature can be adjusted if local velocities are too high 
for the desired thickness. This computed velocity will then serve 
as a guide to the choice of an airfoil velocity distribution, which 
shou3.d be chosen to avoid high velocity peaks and steep negative 
gradients. If the average of the velocities on opposite sides of 
the blade camber line is retained in the modification of the velocity 
distribution computed from the stream-fi lament method, the thickness 
will also be retained. 
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Because of the irrotationality of the fluid motion, the velocity- 
integral or circulation around the airfoil must he equal to that 
around a blade but over a -width equal to one blade space. Therefore, 

y'cp'(s) ds = r = - Vx^2) 

where 

]_ tangential velocity entering cascade 
Vx 2 tangential velocity leaving cascade 
r circulation about airfoil 

This relation places a condition on the assumed velocity distribution. 

If the computations thus far have been made in order to select 
a velocity distribution for the airfoil cascade in a compressible 
fluid flow^ an equivalent velocity distribution for the flow of an 
incompressible . fluid must be determined before the blade shape can 
be computed by any method based on incompressible flow theory. For 
subcritical flows the directions of the incoming and discharge 
velocities are the same for compressible and incompressible flows, 
but for incompressible flow the component normal to cascade axis is 
the same upstream, and downstream. The Karman-Tsien compressibility 
correction (reference 7) ov that of Garrick and Kaplan (reference 8) 
may be applied to the velocity on the blade surface to estimate 
roughly the corresponding Incompressible flow velocity distribution. 
The resulting velocity distribution in any case must satisfy the 
circulation condition. This procedure does not give an exact solu- 
tion for compressible flows, but the resultant compressible flow will 
have approximately the desired characteristics of low pressure 
gradients and no high velocity peaks. 


Computation of Airfoil Shape from the 
Chosen Velocity Distribution 

The numerical computation of the quantities involved in the 
preceding analysis, particularly the function is extremely . 

laborious when tables of f2(3, So) are used. Most of the comx)utations 
are therefore executed graphically. In the cascade example, it was 
assumed that the air entered the cascade at an angle of 45*^ from, the 
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cascade axis and was discharged at an angle of -30° from, the cascade 
axis. The assumed velocity distribution is given in figure 3(a). 

The value of the lift coefficient for this airfoil is 3.1. The 
shapes of the Isolated airfoil and the airfoil in cascade are 
computed by the following steps; 

1. Curves for constant f^_' for the isolated airfoil, or 
constant fg (fig. l) for the airfoil in cascade, are drawn. This 
diagram should be made on some transparent material that will neither 
change size nor shape. The coordinates of the curves for constant fg 
are given in table I. 

2. A desired velocity <P' (s) is chosen as a function of the 
arc length of the airfoil (fig. 3(a)). An airfoil shape having the 
desired total arc length is assumed and is drawn to the same scale 
as the plot of f]_ or fg. The drawing is made on grid paper and, 
in the case of a cascade, the x-axis is taken along the cascade 
axis (fig. 4). 


3. The velocity distribution cp'(s) is integrated to obtain 
the velocity potential cp(s). This function is plotted on the same 
chart as the assumed airfoil shape for the corresponding y-coordinate, 
as shoTO in figure 4, by plotting both cp and the y-coordinate of 
the airfoil against s on a supplementary graph. In regions of the 
airfoil where y varies little with s, that is, where the airfoil 
boundary is parallel to the x-direction, cp should be plotted 
against x in the same manner, as shown in figure 4, 


4. In order to find the stream, function at a point ( x, y) 

on the airfoil, a plot of f 2 (s, Sq) aa a function of c-p(Sq) must 


be obtained to evaluate the quantity 


f 2 (s, Sq) dcp(so) of equa 


tion (8). If the chart of f 2 is superimposed on the airfoil with 
one vortex center overlaying the point (x, y), the value of f 2 
may be read at (xq, y^) and the corresponding value of (?>(Xq, y^) 
may also be read from the plot of cp(xq, y^). The value of f 2 (s, Sq) 
is the same as would have been obtained by centering the chart on 
(^o, Yq) because of the symmetry of the function, A succession of 
values of cp and f 2 are obtained in this fashion for various posi- 
tions (Xq, y^) that intersect the f 2 contours, and a plot of 
these points (f 2 , cp) may be made for the assumed position (x, y) . 
This procedure is illustrated in figure 5 for a particular point 
(x, y) on which the f 2 chart is centered. The readings for a 
particular (Xq, yp) are shown by the arrowed lines. The 
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points 1 to 6 on the hlade are shown on the corresponding f 2 curve. 
The discontinuity of between points 1 and 6 is the circulation. 
The discontinuity between 4 and 5 represents the region where fo 
approaches 

5. The proper method of integration then proceeds from 1 through 
6 to 7 and then to the origin, with constant f 2 from 4 to 5. The 
region from 4 to 5 with fg approaching - t» is computed by equa- 
tion (6) or (6(b)); the constant a is assumed to be the radius of 
the near-circle, which corresponds to the value of fg where the 
discontinuity from 4 to 5 occurs. 

The total area including this small addition is 

J cp'(so) f 2 (s, Sq) ds =Y'f 2 (s, Sq) dcp 

which is the stream function due to vortices on the entire set of 
airfoils in cascade. • Where fg = 0 at the points A, B, C, and D 
(fig. 5), the values of cp are noted as fPy^(s), (p^(s), cpQ(a), 
and cpjj(s). These values are used in computing the stream-function 
change caused by rotating the blade. The stream function at the 
point (x, y) may now be computed from equation (8) or (5), and 


. -Vyc + V^y 

A plot of the stream function (variation from the mean value) 
is shown in figure 6 for the initially assumed shape. Corresponding 
points on adjacent airfoils have a difference of ini^/Yyd. equal to 1.0. 

6. When \l» (a) is known at a sufficient niunber of points, the 
airfoil may be rotated as previously described. For the isolated air- 
foil, equations (14) and (17), may be used directly. For the airfoil 
in cascade, the coordinates of the centroid of the airfoil must first 
be computed by 


2:c 



cp'(so) dSo 
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Before equation (17) can be used to compute d\|//dp, the variable 
quantity ^x.r be computed. Tlie vortices in the region fg^ 0 

are considered to be uniformly distributed along the cascade axis 
and the velocity induced by such a distribution, is 



where 7 is the vortex strength per unit length along the cascade 
axis for fg> 0. Therefore, 


1 ^ 

^x = m ^3o 

where the integral is to be taken over the regions fg> 0. The 
region fp>0, Jq> J contributes . a positive component to v^, 
whereas the region fg>0, Yq<j contributes a negative component. 
The computation Is simply carried out by making use of the fact that 
the integral for v^ is the difference between values of cp at 
points where fg = 0. The values of cp^(Sq), q?g(sQ), c;^(Sq), 
and cPj)(Sq) from step (5) are used at this point to obtain 


2 v^d =1 cp'(Sq) dSQ = cp^ - cpjj + r - (cp^ - c?^) (18) 

u 

where r is introduced because o f the discontin uity in cp at the 

trailing edge. The term cp^r^ “ ^ gives the effect of the 

vorticity in the region fg(s> near the trailing edge, and 

the term. cpQ -cpg gives the effect of the vorticity in the 
region fg(s, Sq)> 0 near the leading edge. If either the leading 
edge or the trailing edge lies in the region fg(s, Sq)<0, only 
two points of intersection will remain and one of the two groups of 

terms in equation (I8) will vanish. The quantity / c?'(so) '^®o 

is added to the x-com.ponent of the original uniform stream velocity 
and the quantity d\j//dP of equation (17) may be computed for a 
number of points and the angle p computed from equation (14), using 
the values of (x^, y^) just determined. After these computations 
have been made, the airfoil is rotated through the angle p, and the 

value ^ assigned as the value of the stream function of 

the point after rotation. 
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7, A value of \J/(s) is knovn at points along the airfoil boundary. 
The mean value over the airfoil ^ is subtracted from \1/ leaving A\V. 
For the isolated airfoil, no subtraction is necessary. Each point is 

moved a distance in the direction perpendicular to 

A,/ + V 

the velocity computed in step 6. The curve joining the points in their 
new positions is the adjusted airfoil. 

8. The total arc length of the adjusted airfoil will be different 
from the original one, in general, although local changes in length 
will be negligible. The airfoil is so scaled that the length of the 
suction side is the same length as it had been before distortion 
because this surface is the critical surface of the airfoil. This 
process will result in a change in length of the pressure side. The 
velocity over the pressure aide cp'(s) must then be altered in such 

a manner that the difference in potential between the two stagnation 
points remains the same. As a result, the q[uantities that retain 
specified values are the length and the velocity distribution on the 
suction side and the circulation around the airfoil. The entire 
procedure is repeated with the adjusted shape until the variations 
in the stream function result in very little change in the shape of 
the airfoil . 


Discussion of Examples and Techniques 

For the example being computed, the stream functions obtained 
for the initially assumed shape and the first and seventh approxima- 
tions are plotted against the arc length (fig. 6), which is taken as 
zero at the trailing edge and proceeds counterclocln/ise around the 
airfoil as shown in figure 7. The fact that for the initial: 

shape is positive over the first half of the arc length and negative 
over the second half indicates that it is too thick because the 
required distortion in shape will make it thinner. The change in' 
thickness results in a change in velocity distribution over the pres- 
sure side of the airfoil in order to maintain the desired circulation. 
The velocity that was originally assumed, which is equal to the 
vorticity per unit length distributed on the initial airfoil, is 
shown in figure 3(a) and the velocity over the- final shape in fig- 
ure 3(b). The length of the pressure side has increased and the 
velocity has decreased in the proportion of 1:1.1. 

Over the section of the airfoil that has collapsed to zero thick- 
ness, the surface velocities of figure 3(b) may not have been obtained, 
but the loading (circulation per unit arc length), which is the 
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clii’fe^rence in ulie velocities on opposite sideS; has been realized. 

In practice, this collapse is prevented by , increasing the assumed 
velocity on the airfoil surface. 

If the initially assumed airfoil shape has a thickness that 
differs considerably from, the correct one, the process of shape 
adjustment will converge rather slowly. The labor can be reduced, 
however, by computing the stream function at a few points on the 
airfoil and locating these points to determine the thickness. This 
procedure is followed for the first few approximations until the 
thickness of the airfoil is fairly accurate. The stream function 
is then computed at a larger number of points, particularly near 
the leading edge, in order to get more detail of the shape. 

If a velocity distribution is arbitrarily specified, the 
resulting shape maj'' not be a physicalD.y real airfoil but may result 
in a figure-G or a col.lapsed shajja (zero thickness over a portion 
of the blade). It is then necessary to modify the velocity 
distribution to obtain a real shape; these modifications should be 
selected to keep the desirable properties of the original distribu- 
tion. Velocity peaks and steep velocity gradients, which tend to 
occur on the suction side of an airfoil, are to be avoided. If the 
airfoil collapses, the vorticitles of the two sides tend to cancel 
each other and the remaining vorticity represents the difference in 
velocity across the thin airfoil rather than the velocity along the 
boundary. 

The method was also applied to the design of a thin airfoil 
(camber line) in a cascade. The vortex distribution is eq^uivalent 

-to— load-d-istr-ibution— (-d4-f-ferenGe— i-n— vel-oe-i-ty— across— the— a-i-rfoi- 1 -) 

rather than velocity as in the case of a thick airfoil. The 
velocity diagram for the cascade and the desired velocity difference 
are shown in figure 8. The value of the lift coefficient of the 
resultant airfoil is 4.1. The initial shape was obtained by 
assuming zero spacing between the airfoils. The initial shape and 
the first and third approximations to the airfoil shape are shovm 
in figure 9. There is very little difference between the second 
and third approximations. The third approximation is redrawn in 
this diagram to show the spacing between airfoils. The convergence 
of the method is sho>m graphically in figure 10. The varia- 
tion A\lf of the stream, function from its mean is divided by Vd to 
make it dimensionless and is plotted against the arc length along 
the airfoil where s = 0 at the trailing edge. The stream function 
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computed on the second approxiroation is nearly constant, which gives 
the third approximation almost the same shape as the second one. 

The rapid adjustment of camber contrasts with the slow adjustment 
of thickness. 


Flight Propulsion Research Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, March 4, 1947. 
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APPENDIX A 


TELOCITf DISTRIBUTION ON THE DERIVED AIRFOIL 
AT DIFFERENT FLOW ANGLES 

Confonaal mapping. - When an airfoil is given, the velocity 
distribution over its surface must frequently be found at different 
angles of attack. This problem may be solved by the method of con- 
forcjal mapping, which consists of mapping the region exterior to 
the airfoil on the exterior of a circle. The velocity around the 
airfoil is obtained fixom the known velocity around the circle. 
Procedures for finding the function that maps a given airfoil into 
a circle are presented in references 1 and 9 for the isolated 
airfoil and references 3 and 10 for the airfoil in cascade. 

In general, the procedure for finding the mapping function 
of an airfoil is a laborious one. But when, as in the present 
case, the velocity distribution over the airfoil at a particular 
angle of attack is known, the correspondence between points on 
the airfoil and on the circle, and hence the flow velocity at 
other angles of attack, can be obtained very easily. Indeed, 
the correspondence of points and the velocities for various 
angles of a,ttack can be obtained by the method of V/einig and 
Gebelein (reference 11) from the initial data without knowing 
the airfoil shape, because the complex potentials of the airfoil 
plane and the mapping circle plane are equal. Before the airfoil 
is designed, therefore, it is possible to check whether the air- 
f_oll_t.o_b.e_Gomputed— w4.il— be— 8at-isfaetory--under-conditions~different 
from the design condition. 

Isolated airfoil. - The flow about any airfoil shape can be 
mapped on the flovf about a unit circle in such a way that corre- 
sponding points have the same potential. The flow about the air- 
foil is given and the potential function cp(s) at each point is 
computed. If the potential function on the airfoil is computed 
by integrating the velocity from the stagnation point at the 
trailing edge in a countercloclcwise direction around the airfoil 
oriented like the one in figure 1, the potential will be zero at 
the trailing edge, decrease to a minimum stagnation 

point at the leading edge, and then increase to a value equal to 
the circulation F at the trailing edge. The corresponding flow 
about the circle is determined by the conditions that it must 
have the same values of ^ ^ correspondence to 
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exist between all airfoil and circle points. If central 

angle of the stagnation point on the circle that corresponds to the 
trailing edge of the airfoil, 

= —(cot + 0Q1 + n/2) (Al) 

Equation (Al) can be solved numerically for 0|p because all the 
other quantities are taiown. The velocity at infinity in the circle 
plane Vq can then be determined from the Kutta-Joukowsky condition, 
which requires that 0^ be a stagnation point; that is, 


4jt sin 0qi 

The velocity potential at points on the circle is 


Cp ss •2V- cos 0 + rr— + 
c c 2n u , . 


27 T 


(a2) 


(A3) 


r> 

The quantity 2Vq cos " 2n ^ constant that is subtracted 
in order to make cp^ = 0 at the stagnation point corresponding to 
the trailing edge . 


The corresnondence of points on the airfoil with points on the 
circle is obtained by associating points where cp(s) = cp^. The 
velocity on the circle at a uniform stream flow angle is 

Vc = 2Vc I 8in(0 + a]_) - sin (0 qi h- cxq^)J (A4) 

The nature of the conformal transformation is such that the ratio 
of the velocity at a point on the airfoil to the velocity at the 
corresponding point on the circle is independent of angle of attack. 
Therefore, the -velocity cpj_'(s) on the airfoil at flow angle is 


cpi'(s) cp'(s) 


(AS) 


Vc(0,Oi) v^,(0,O) 
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where the design flow angle is taken as zero. Equation (,A5) can 
be used to compute the velocity distribution on the airfoil except 
at the tvro points that v/ere stagnation points at the design angle 
of attack. 

Airfoils in cascade. - The flow about a cascade of airfoils 
can be mapped conformally into the flow about a unit circle with 
two singular points located on the real axis symmetrically with 
respect to the center of the circle. These singular "points cor- 
respond to the points at infini'by in front of and behind the 
cascade, respectively. In a cascade of airfoils, the distance 
of these points from the center of the circle is uniquely deter- 
mined by the same conditions that determine the flow about the 
circle in the isolated case; namely, the cii'culation per airfoil, 
the velocity "potential at the leading edge, the blade spacing, 
and the upstream and downstream f lovr angles . 

The distance from the singular points to the center of the 
circle is denoted by e^. The flow about the circle is such 
that the location of the stagnation -points 0qi is determined 
by the relation 

r sin 0m cos 0m 

QQS X -f sin A (A6) 

2"Vd sinh K cosh K 

vjhere X is the angle of inclination of the mean stream to the 
normal to the cascade axis. . (See reference 6 for details). The 
quantities P, "V, d, and X are known from the flow in the 
cascade plane and therefore equation (A6) provides a relation 
-■"dDetween" “K"~an"d“t"he~ln'^c"ation"of ~the" staghat'ioh^podnts’^ 

The velocity jjotential at any point on the circle is 


cpQ Q = ~i sin X tan"^ - cos X tanli"^ - 

n \ sinli K cosh K ) 

r ^ _i tan 0 i Vd (' . ^ ®T ^ , -1 "^T^ 

— tan j — [ sin ?\ "can cos X tanh ) 

2rt tanh K L it \ sinii K cosh K/ 


r tan 0qi 

+ — tan”""^ , 

2it tanh K J 


(A7) 
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The expression in brackets is a constant so chosen that the potential 
will vanish at the stagnation point corresponding to the trailing 
edge of the airfoil. In order to map the cascade on the circle, it 
is required to find K so that the value of cp^ ^ at the 

stagnation point corresponding to the leading edge of the airfoil, 
is equal to cPjnj[n> value of the velocity potential there. In 

s'uccessive aipproximations, a value of K is assumed and equation (A6) 
is solved simultaneously with the identity 


/ sin 0 r 9 I cos 0 P 

viiSTKy + v^ihicy ^ 


for the variables 
equations are 

r 


sin 0 
sinh K 


and . 


cos 0 
cosh K* 


The solutions of these 


sin 0g “2^ cosh^K sin X±cos X>y4o3h^K- cos^X cosh^K sinh^K 


sinh K 


cosh'^ K - cos^ X 


cos 0g 

and obtained from equation (A6) 

cosh K 


/sin 0JJ cos 0 




N 


) 


and ( - 


sin 


Qn 


COS 




sinh K" cosh K 
tion XA^O) 'to find the value of 
not equal to cPminJ 


. These solutions, 
are substituted in equa- 
at 0 = 0i\T. If 


^sinh K cosh 

^c,c 

another value of 

that cpg decrease as K 

is evaluated, care should be taken to use consistent values of the 
inverse tangents. After two values of K and cpg c(%) deter- 


^c,c(%) 

K is chosen, on the premise 
is decreased. \<dien 


mined. Interpolation or extrapolation may be used for new values 
of K. 


When IC has been found, it is used in equation (^Af)) to eval- 
uate cpg Q at values of 0 all around the circle. A point on the 
circle corresponds to the point on the airfoil where cp(s) = 

The velocity at the point 0 on the circle is 


c , c 


(0,a) = H 


sinh 2K 


It cosh 2K - cos 20 L 


I cos X 


cos 0 
cosh K 


cos 0ip\ 
cosh kJ 


/ sin 0 sin 0^ \ 
+ sin ainh K / 


-4A0) 
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and the velocity on the airfoil at any other mean flow angle 
ia 


as in the case of the isolated airfoil. 


(A-9) 


The designed airfoil was mapped on the unit circle hy the 
method described. The constant K, the natural logarithm of the 
distance from the singular points to the center of the unit circle, 
is 0.075. The correspondence of points on the airfoil with those 
on the circle is plotted in figure 11, which shows the arc length 
of the airfoil as a function of the central angle of the circle. 

The velocity at any point on the airfoil for any angle of attack cu 
may be obtained from equations and (A9) , using the velocity 

distribution as in figure 3(b) and the relation between s and 9 

as in figure 11. The ratio is equal to || (radians) 

and need be coraputed only once for any given airfoi.l. 
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APPSHDIX B 


DERIVATION OF THE CASCADE EQUATION 

An eq,uation is to be developed for the complex velocity at any 
point in the field of flow of a fluid past a row of equally spaced, 
congruent bodies. Coordinate axes are chosen with the origin inside 
one of the bodies and the x-axia in the direction of the row. 

(See fig. 12.) The body containing the origin is denoted by Bq^ 
bodies along the positive direction of the x-axis by B-j_, B2, etc., 

and along the negative direction of the x-axis by B_q_, B_2^ etc. • 

A circle A of small radius is drawn about the point z where the 
velocity is to be determined. A rectangle E is drawn with its 
center at the origin and its sides parallel to the axes of length 
(2N+l)d and width 2t, which contains- the bodies B_jj . . . B_p, 
Bq, B]_ . . . B|^, and the circle A. If a side of the rectangle 
intersects one of the bodies, the side may be distorted to go 
around the body with no essential change in the proof. The func- 
tion w (Zq)/zq-z is an analytical function of z in the i*egion 
inside the rectangle R but outside the bodies Bn and the 
c ircle A . 


Therefore 



N 


\ 




I ■; — 


J "0 


= 0 ^:^Bi) 


n = -N B. 


n 


The first integral can be broken up into four integrals, one 
along each side of the rectangle, namely. 


r 


pXN+l/2)d 


Zq-Z " 


rt 


^y-(N+i/2)d ° 

-(N+l/2)d 


w'(Xo-it) 

dx_ + 

X_-lt-Z O 


t;-t 


w'[(N+l/2)d+iy^]_ 
“(Nri72)d+iyo'^ ^*^^0 


A“t 


(N+l/2)d . 


w'(x^+it) 
Xg+it-z ^o 


w’ [-(N+l/2)d.+iyQ] 

-(N+172 )d+iyQ-~ ^^^o 
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In an evaluation of these integrals, the function w'(zo) is 
periodic, with period d, and approaches a constant value 
infinitely far from the cascade; that is. 


w’(Xo + iyo)~^W2’ 

and 

w'(xo + iyo)— 

From the last of these conditions, it 
w'(Xq -it) = W3* (x^ 


as -^ co 

as yo — > -» 
follows that 
* it) + w^' 


where 


Wj' (x^ - it) as t -7» 


Therefore, the first integral on the right side of equation (B20) 
is 


,''^(N+l/2)d 

/ w'(xQ-it) 

I Xo-it-z 

J-(N+l/2)d 


wi 


/<N+l/2)d 


J-(N+l/2)d 


Xp-it-z 


+ 


,->(N+l/2)d 

I 

» 

I 

f 

J-(N+l/2)d 


W3'.(xQ-it) 

Xg-it-Z 


dx^ 



The first of these integrals is 


w 


1 


rv(N+l/2)d 

i 


,j-(N+l/2)d 


Wi 


log- 


[(N+l/2)d-it-z] 

[-(N+l/2)d-it-z] 


— ^ rt iwj^ ' 


as N — )<» and t— 4 oo^. provided that ^o. The last integral 
in equation (B3j)is 
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/.(N+l/2)d 

/ W5’(xQ-it) 

■ ^ Xo-it-z 

J -(W+l/2)d 


dx. 


N_ j'(n+l/2)d. 

^ " ! W5'(xo-it) 

n = -N 'J (n-l/2)d 


dx. 


\" 
= z 


n = -W ij 


ni/2d 

~l/2d 


W3'(Xo~lt) 

x^+nd-it-z 


nl/2d 

W /^l/2d 

/ ^3'(^o-^~^) 


1 Xo-it-z 

'• -l/2d 

+ y j 

n = 1 vi -l/2d 


2(xo-it-z)w3' (xQ-it)dXQ 
(xQ-,it-z)2 -ri2 d2 


If t is chosen sufficiently large so that |w„‘(x -it)|< e, where 

is any preassigned positive nurober, the integrals are less than or 
equal to 


>l/2d 


rl/2d 


dx^ 


2 XQ-it-z } dxQ 


' - Z.,.y / , i(x -it-z)2 -n2 d2j 

-l/2d ' n = 1 J -l/2d ' ° ' 


<€ 


■'l/2d 


dx 


w .^ ra/2d 

\ I 2 V(xo-x)2+(yo-t-y)2dXc 
+ / 


'\/(x-x)^+(y I (x -x)2+(y -t-y)2-n2d2 


[J -l/2d " n=li’-l/2d ° 

When W-~^®^, this quantity aTJproaches 


rl/2d 

e I I cot J- + (yo-t-y)2 dx^ 

/ vL Q. 

J -l/2d 
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This integral is finite and, hecause e can he made arbitrarily 
small as t — > op ^ the last integral in equation (B3) approaches 
zero . Therefore , 


(N+1/2 )d 

/ 

/ 


r 

J. 


(N+l/2)d 


w'(xQ-it) 



itiWj^' 


as N s-®* and t rco wnile _ —^0. in the same way and under 

N 

the same conditions, 


^ N+1/2 )d 

/ 

i 

I 

(N+1/2 )d 


w’(x^+it) 

x^+it-z 



The second and fourth integrals on the right side of equa- 
tion ,'0B2) can be evaluated by combining them. Because w' is 
periodic, 


w' [ (N+1/2 )d+iyQ] = w' [ -(N+1/2 )d+iyQ] 


and therefore. 



w’ C(N+l/2)d+iy^^]^^ 
(N+1/2 )d+iy^-z 


w'[ -(N+1/2 )d+iyQ] ^ 
- (N+1/2 )d+iyQ-z 

'./ 1 


i 



-2(N+l/2)d w' [(N+1/2 )d+iyo] 
(iy^-z)2 - (N+1/2)2 d^ 


idyo 


The velocity w' [(N+l/2)d+iyo] is bounded for all values of 
that is, there is a constant W such that Iw' [ (N+l/2)d+iyQ]j<W. 
The absolute value of the integral is less than 
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2d(N+l/2) W 




(iyQ-z)2 - (N+1/2)2 ^2 


"S 2d(N+l/2) W 


Ai 


-t 




+ (N+1/2)2 d2-x2 


_ 2d(N+l/2) W 


.V&+1/2 ) 


2 d2-x2 


tan"^ 


t-y 


VSi^/2) 


— tan 


-1 :t-y 


2 ^ 2.^2 


/\/(N+1/2)2 d2-x2 


As N — >•<» and t — >oo^ this quantity approaches zero. It has been 

shown, therefore, that when t — and — >0, 

" ' N ' . 


E 


v’(zq) 

Zq-Z 


dz. 


«i (w^’h- w^') 




By the residue theorem, 



2jt iw ' ( z ) 


:t«5) 


The periodicity of w'(z) implies that 
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n = -N v/ 
B 


^ dz^ 

2q-z o 


/ , / z +nd-z o 

n ° 


r X? r ' r 

/ ^’(zq) . \”' h^’(Zo) 2(zo-2) , « s , 

' Zq-z J (z -z)2 -„2 42 ° ^ d ° d ° ° 


n = 1 B, 


,i(B 6 ) 


as N — 


When eq,uations .(B4)^ (fB-S), and are substituted into 

equation (Bl), the expression for the complex velocity is obtained: 


w'(z) = |(vi'+W 2 ') “271 j 7 ^' (Zo) cot I (zp-z) dz, 

®0 


4<B7) 


The conjplex potential is obtained from eq_uation .(®7)) "^7 integrating 
with respect to z and neglecting the arbitrary constant, 


w(z) = z / w’ (z^) log sin J (z-z^) dz^ (]fB 8 ) 


Wi * + Wg ^ 

where w^q* = is the mean stream velocity. 

Ct 
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APPEWDIX C 
SYMBOLS 

The main symbols used throughout the report are listed here 
for convenience of reference. 

d distance between successive airfoils in cascade 

fl= ^ Ux-Xo)2 + (y-yo)2j 

f 2 - ^ i sin^ 1 (x-Xq) + sinh^ | (7-yo)J 

K natural logarithm of distance fi-om singular point to 
center of circle corresponding to cascade airfoils 

I total arc length of ^airfoil 

s arc -length parameter corresponding to z 

Sq arc -length parameter correspondong to Zq 

Vq local velocity on circle corresponding to isolated airfoil 

^c,c Local velocity on circle corresponding to airfoil in cascade 

v^ velocity induced by vortices in region f 2 > 0 

V magnitude of uniform or mean stream velocity in airfoil or 

cascade plane (fig. 1) 

Vjj magnitude of uniform stream velocity in circle plane 
x-component of uniform or mean stream velocity V 
V„ T. resultant local mean stream x-component of velocity V 

A, X 

Yy y-component of unifoim or mean stream velocity V 
w complex potential function, (cp+ i'l') 

w ' Comdex velocity of uniform stream for isolated airfoil, 
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Wjjj' coE^plex velocity of mean stream for airfoil in cascade 
‘>m‘ = |-('^l’ + W2') = - iVy) . 

^c-»^c coordinates of point about vhich airfoil is rotated (centroid 
of vortex distribution for cascade airfoils) 

2 coominate of ijoint where stream function is computed, 

(x + iy) 

Zq coordinate of point where vortex is located, (x^ + iy^) 

<x angle of inclination of uniform stream velocity to x-axis 

3 angle through which airfoil is rotated 

7 (Zq) vortex strength per unit arc length at Zq 

r circulation about airfoil 

6 central angle of circle 

0JJ angle of stagnation point on circle corresponding to leading 

edge of airfoil 

0J, angle of stagnation point on circle corresponding to trailing 

edge of airfoil 

A angle of inclination of mean flow to normal to cascade axis 

(fig. 1) 

qj velocity potential on airfoil, E [w(z)] 

cp^ velocity potential on circle corresponding to isolated airfoil 

% c velocity potential on circle corresponding to airfoil in 
’ cascade 

fp^ values of cp at points A, B, C, D, respectively, where 

cpQ^rpp^ the curve of intersects f 2 (s,SQ) = 0 (See fig- 

’ ure 5.) 


'■?min 


velocity potential at the leading edge of the airfoil 
stream function, I[w(z)] 
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•\l/^ stream function of mean stream of cascade flow 

stream function of uniform stream flowing about isolated 
airfoil 

mean value of stream function over airfoil 

variation of stream function, - 'I' ) 

Subscripts 1 and 2 when appended to w', V, and indicate entrance 
and discharge values, respectively. 
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table I. coordinates or 


(a) Values of — 
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0.05 

0. 10 

0. 15 

0. 20 

0.25 

0.30 


.35 

0.40 

0.45 

0. 50 


.048 I 


k .28 I .O545I0.O229 


.0618 
I- .24 I .0699 
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Entrance velocity diagram 


Reflected velocity diagram 


Figure I. — Notation for cascade flow. 





NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS 



Figure 3. - Prescribed velocity distribution for thick airfoil in cascade 



Fig. 4 
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Figure 4, - Plot of airfoil and velocity potential for use In computation. 
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Fig. 5 
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Figure 5. - Superposition of figures 2 and 4 to obtain plot of fj against <e. 



Figure 6. - V 4 riation in stream function along initial shape and first and seventh approximations of airfoil cascade 
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Fig. II 



Tioure II.:- Correspondence between points on airfoil and points on unit circle 

by conformal transformation. 
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Fig. 12 
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Figure 12, - Diagram for derivation of equation of flow about cascade. 



